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Abstract. We apply the discrete multiscale expansion to the Lax pair and to the first few 
symmetries of the lattice potential Korteweg-de Vries equation. From these calculations 
we show that, like the lowest order secularity conditions give a nonlinear Schrodinger 
equation, the Lax pair gives at the same order the Zakharov and Shabat spectral problem 
and the symmetries the hierarchy of point and generalized symmetries of the nonlinear 
Schrodinger equation. 

1. Introduction 

Reductive perturbation techniques [19,20] have proved to be important tools for finding 
approximate solutions of many physical problems, by reducing a given nonlinear partial 
differential equation to a simpler equation, often integrable [3], and for proving integra- 
bility [3-5,10,21]. Recently, after various attempts to carry over this approach to partial 
difference equations [1,11,13] we have presented a procedure for carrying out a multiscale 
expansion on the lattice [7,12,14] which seems to preserve the integrability properties [8]. 
To get a better understanding of the application of the reductive perturbation technique 
on difference equations, after an introduction in Section [2] on multiscale expansions on 
the lattice potential KdV equation (lpKdV), we discuss in Section [3] its application to the 
spectral operator, as was done by Zakharov and Kuznetsov in their pioneering work in 
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1986 [21] for the KdV equation. Later on we apply, in Section HJ the multiscale expansion 
to the symmetries of the lpKdV [15]. Section is devoted to a few conclusive remarks. 

2. Multiscale expansion on the lattice 

The aim of this Section is to give a terse survey on the multiscale analysis on the lattice 
and its application to the reduction of the lpKdV. We refer to [7, 12, 14] for further details. 

2.1. Shift operators defined on the lattice. Let u n : Z — > R be a function defined on 
a lattice of index ?i6Z. One can always extend it to a function u(x) : R — > R by defining 
a real continuous variable x = na x , where o~ x G R is the constant lattice spacing. 

An equation defined on the lattice is a functional relation between the function u n 
and its shifted values w„±i, u n ± 2 , etc, expressed in terms of a shift operator T n such 
that T n u n = u n+ i. 

For the continuous function u(x) we can introduce an operator T x , such that T x u(x) = 
u(x + a x ). The Taylor expansion of u(x + a x ) centered in x reads 

OO j 

T x u{x) = Y,^{x), (1) 

where m«(i) = d 1 u{x)/dx l = d l x u(x), with d x the total derivative. Eq. ^TJ) suggests the 
following formal expansion for the differential operator T x : 

rp _ cr x d x _ \ s u x M 

J-x - e - i\ a v 

i=0 

Introducing a formal derivative with respect to the index n, say 5 n , we can define, by 
analogy with T x , the operator T n as 

T n = e s - = T-T- ( 2 ) 
The formal expansion can be inverted, yielding 

00 f_ iV~i 

S n = lnT n = ln(l + A n ) = V r — A^, (3) 

z — ' 2 

t=l 

where A n = T n — 1 is the discrete right difference operator w.r.t. the variable n (i.e. A n u n = 

u n+l ~ u n)- 

Following [12, 14] we say that slow-varying function of order I iff A^ +1 u n = 0. 

Hence the 5 n operators are formal series containing infinite powers of A n , but, acting on 
slow-varying functions of order £, they reduce to polynomials in A n of order at most t. 

2.2. Dilations on the lattice. Let us introduce a second lattice, obtained from the first 
by a dilation. For x G R we can visualize the problem as a change of variable between 
x and x\ — ex, < e <C 1. On the lattice this corresponds to a change from the index 
n = xja x to the new index n\ = X\/o~ Xl , where a Xl is the new lattice spacing. Assuming 
that cr Xl ^> cr x we can set a x = ea Xl , < e <C 1, so that n\ = een. As n, n\ G Z, ee is 
a rational number and one can define in all generality ee = M±/N 1 with Mi, JVel 
However, if we want the lattice of index n\ to be a sublattice of the lattice of index n, we 
have also to require that M x /N = 1/M with Mel 



MULTISCALE EXPANSION AND INTEGRABILITY PROPERTIES OF THE LPKDV 



3 



The relation between the discrete derivatives defined in the two lattices is given by 
[7,9,13,14] 

oo 

AX = j!^^A> ni . (4) 

1=] 



The coefficients Py read 



Mi x 



fc! 



where iSf and are the Stirling numbers of the first and second kind respectively. 

If u n is a function of infinite order of slow-varyness, i.e. i = oo, then Eq. (j3J) implies that 
a finite difference in the discrete variable n depends on an infinite number of differences 
on the variable n\. 

2.3. Discrete multiscale expansion. Let us now consider u n = u n]ni as a function 
depending on a fast index n and a slow index n\ = n(M\/N). At the continuous level, 
the total derivative d x acting on functions u(x; x\) is the sum of partial derivatives, i.e. 
d x = d x + ed Xl . As 

we can write the total shift operator T n as 

j n _ e 5„ e (Mi/AT)a ni _ T n T( Ml/N \ (6) 

where the partial shift operators T n ,T ni , defined by T n u n;ni = u n+ i ;ni and T ni u n - ni = 
u nm+1 , are given by 

T — V" ^ T^'h/ N ) - V s ( M i/ N Y h 

1=0 ' i=0 

and <5 ni is given by Eq. (j3J) with n substituted by nj. 

Eq. (0) can be extended to the case of K slow variables Xj = e 4 x, 1 < i < if . Then 
the action of the shift operator T n on a function u n .{ n .\K_^ depending on both fast and slow 
variables can be written in terms of the partial shifts T n , T Hi as 

K 

T n = T n \\rtr\ (7) 

where the e ni 's are suitable functions of e and e depending parametrically on some integer 
coefficients G N, 1 < i < K. 

To carry out the multiscale expansion of the fields appearing in partial difference equa- 
tions with two independent discrete variables, one has to consider the action of the oper- 
ator ([7j) on a function depending on two fast indices n and m, and on a set of K n + K m 
slow variables {ni\ 4 J\ and {mAf^, (we shall use the notation u s x k„ , X K m for such 

L Jt-L *- J l -L \ 71,771, \Tl.i j |_ ^ , \TH,i j ^_ ^ 

functions). Notice that in principle it is possible to consider K n = K m = oo. We assume 
a common definition of the small parameter e for both discrete variables n and m but we 
denote with Mj the integers for the slow variables n« and with Mj the ones for m,. We 
have: 

Mi M 
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Hereafter we shall assume K n = 1 and K m = K. 

2.4. Multiscale expansion of the lattice potential KdV equation. The lpKdV is 

given by [17]: 

[fjj(T n T m 1) -\- (^(T n 7m)]fii,m (-^n T m ^)u ntm (T n T m l)tt njm 0, (8) 

where \i = p — q and ( = p + q, and p,q, p ^ q, are two real parameters. The linear part 
of Eq. ([8]) has a travelling wave solution of the form u n ^ m = exp {i[nn — u(K,)m]} with 

( C + I 1 K \ 

u>(k) = —2 arctan tan — I . (9) 

According to [7] the multiscale expansion of Eq. ([8]) is performed taking into account 
that 

oo 1 

un, m = E £ w^ ){ni > wf=i)e ia( ™\ = (io) 

aeZ k=l 

The following statement, proved in [7], provides the multiscale expansion of the lpKdV 
([8]) at the lowest orders of 1/N. 

Theorem 1. The multiscale expansion of Eq. (EP gives the following results: 

(1) 0{l/N): 

• o = 0: the equation is identically satisfied. 

• a = 1: one gets a linear equation identically satisfied by taking into account 
the dispersion relation 

• \cx\ > 2: one gets a linear equation whose only solution is u ™ = 0. 

(2) 0{l/N 2 ): 

• a — 1: one gets a linear equation whose solution is 

uf> = «i (n 2) {mi}^), n 2 = n 1 Tmi, (11) 
provided that 

/-2 2 

M x = TS{fi- Ce iK ) , Mi = Se iK ±-^-. 

v ' \ie 1K — C 

Here S = rexp(i9), with r > and 6 = — arctan [(£ sin k)/(C cos k — n)}, 
assures that Mi and Mi are positive integers. 

• o = 0: one gets 



Se*«(ji + t) (/i-Ce L 

where uf^ = Ui(n 2 , {mi}f =2 )- 
a = 2: one gets 

(2) , (i)x 2 1 + e 1K 

u\ = t 2 {u\ ') , r 2 



(3) Oil/ TV 3 ); 



;i-e-)(/i + C) 
where u 2 = u 2 (n 2 , { m i]f= 2 )- 



MULTISCALE EXPANSION AND INTEGRABILITY PROPERTIES OF THE LPKDV 



5 



• a — 1: one gets the following (defocusing) dNLS: 

iS m2 u? ^ Pl S 2 n2 u? + p 2 u?\u?\ 2 , (12) 

where 

p( r<2 (.( 2 — P 2 ) s i n K %(p(( — p) (1 + cos K ) 2 sm K 

Pi = = , P2 = — o- 

M 2 (C 2 + p 2 -2(p cos k) M 2 (p + ()(( 2 + p 2 -2(pcosk) 2 

• a = 0: one gefe 

S n2 u ( 2 0) = n (u^u^ + u^u^J - r 3 (u^5 n2 u^ - u^8 n2 u^A , 
with 

2i sin k 

/i + C 

w/iere u^ 0) = 4 0) (n 2 , {mi}f =2 ) and u[ l) = u£\n 2l {m,}£ 2 ). 

• a = 2: one gets 



uf =T i U ( t\5 n2 U^ ) )+2T 2 U ( S2 ) ; 

where = uf\n 2 , {mi]f =2 ) . 



2Se iK (a + l3e 1K ) 

T4 ~ ( e -- mp+cy 



We have given above just those results necessary to get a discrete nonlinear Schrodinger 
equation (dNLS) as a secularity condition and its symmetries. 

3. MULTISCALE EXPANSION OF THE LPKdV SPECTRAL PROBLEM 

As shown in [14] there are many forms for the linear problems associated with the lpKdV. 
The first to be introduced [17] is given by first order 2x2 matrix difference equations. 
Later on [14] it was shown that the matrix Lax pair could be easily reduced to a scalar 
non-symmetric difference equation of second order, used by Boiti et. al. [2] to integrate an 
alternative form of the equations of the Volterra hierarchy. In [15] it was moreover shown 
that by a Miura transformation it is possible to associate the lpKdV with the Toda spectral 
problem introduced by Manakov and Flaschka [6] when the field b n (t) = 0. 

One could start from any of the three linear problems defined in the previous paragraph 
to do the multiscale expansion. However we choose as starting spectral problem the one 
whose second derivative is expressed in a symmetric form, i.e. the discrete Schrodinger 
spectral problem used to integrate the Toda and Volterra equations. 

The n-evolution equation of the (scalar) spectral problem of the lpKdV ([H]) may be 
written as [15]: 

0n-i + a n (f) n+ i = n(j) n , (13) 

with 

Ap 2 

an = [2p - (T 2 + lK, m ] [2p - {T n + T-> n , m ] ' 

Here p G C is the spectral parameter. 

Our aim is now to perform the multiscale expansion of Eq. ( jTBl) in order to get the 
corresponding evolution equation of the spectral problem of the dNLS ( jT2l . We refer 
to [21] for the continuous counterpart of this analysis. 
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To expand Eq. ({TBI we consider the development ({TO]) for the field w n>m , with the re- 
striction (II ip . while the function <p n will be expanded according to the formula: 

OO 

^« = EE w^ ){n2 ' wf= 2 ) eiQ( ™ )/2 > = (14) 

a odd k=Q 

At order Oil), the multiscale analysis of Eq. f|T3|) suggests the following expansion for 
the spectral parameter \i: 

OO 

„ = 2cc. (15) 

Taking into account Eq. ( 1T51) we proceed to the order 1/N of the multiscale expansion 
of Eq. flJJ). We have: 

o (!) 

(1) , 2U 1 2 ( K \ 7(1) _ ^1 ,(1) 



4*^ + ^~ COS* ^ ^ = (16) 

for a = 1. The corresponding equation for a = —1 is given by performing the complex 
conjugation of Eq. ({TBI) . The coefficients of the higher harmonics in Eq. (fT4|) can be written 
in terms of (j) G . For instance, for a = 3, we have: 

(3) _ e +e (1) (1) 

<Pi - 1 _ giK % </>0 > 

By a proper rescaling of 0q and fi\ Eq. ( {TBI) is equivalent to the standard Zakharov- 
Shabat spectral problem of the integrable NLS [18]. 

4. Multiscale expansion of the first two generalized symmetries 

Lie symmetries of a lattice equation H>(u ntTn , T^u n ^ mi T^u ntm , . . .) = are given by those 
continuous transformations which leave the equation invariant. From the infinitesimal 
point of view they are obtained by requiring the infinitesimal invariant condition 



pr l„ jm D 

where 



0, (17) 



By prX n/m we mean the prolongation of the infinitesimal generator X n ^ m to all points 
appearing in D = 0. 

If F n m = -F n m (-u„ im ) then we get point symmetries and the procedure to get them from 
Eq. (TTTT) is purely algorithmic [16]. Generalized symmetries are obtained when F n m = 
T^u ntm} . . .). In the case of nonlinear discrete equations, the Lie point 
symmetries are not very common, but, if the equation is integrable and there exists a Lax 
pair, it is possible to construct an infinite family of generalized symmetries. 

In correspondence with the infinitesimal generator ( {TBI) we can in principle construct a 
group transformation by integrating the initial boundary problem 

dUn,m( X ) _ r /- /U r r +~ ,>x r r ±~ 



K,m(Un,m(^),T n U n , m (X),T m U njm (X),...), M„, m (A = 0) = U n , m , (19) 
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where A G K is the continuous Lie group parameter. This can be done effectively only in 
the case of point symmetries, as in the generalized case we have a differential- difference 
equation for which we cannot find the solution for a generic initial data, but, at most, 
we can find some particular solutions. Eq. (iT7j) is equivalent to the request that the A- 
derivative of the equation D = 0, written for £t nj?n (A), is identically satisfied when the 
A-evolution of £t n>m (A) is given by Eq. (TT9l) . This is also equivalent to say that the flows 
(in the group parameter space) given by Eq. ( fT9"l) are compatible or commute with D = 0. 

In [15] one can find an infinite hierarchy of integrable generalized symmetries for the 
lpKdV (jHJ) constructed by looking at the isospectral deformations of the Lax pair. The 
first two symmetries of this hierarchy are given by 

du n ^ m I I (20) 



dX 2p + (T- - T n ) 



dX [2p + (T~ - T n )u njZ 



1 

+ 



4?- < 21 > 



_2p + (1 - TZ)u„ im 2p + (T~ 2 - 1)S„ 

The constant terms appearing in the r.h.s. of Eqs. (I20ll2ip ensure that the above flows go 
asymptotically to zero as u n , m — > cost. 

To perform the multiscale expansion of the generalized symmetries (I20|l2ip we consider 
the following development for the field u n>m , see Eq. f lTUj) : 

<m = EE ]^4 Q) K {m}? =2 , {\i}t*V a{K ™\ 4~ a) = (22) 

a&L k=l 

where Aj = X/N l are the slow- varying group parameters, n<i is given by Eq. fill I) and 

^n,m({Aj 0}j_g) U n ^ m . 

Since Eq. (1T21) involves the harmonic we are actually interested just in those equa- 
tions, arising from the multiscale expansions of the symmetries ( 1201121]) . which are written 
in terms this harmonic. The following statement holds. 

Theorem 2. The multiscale expansion up to order 1/N 4 of the symmetry ( Tj?7j) gives the 
following symmetries for the dNLS [Tfy) ( after a reparametrization of the group parame- 
ters): 

0(1/N) : ^L. = ifi « (23) 
0{l/N-): 2^- = S ^\ (24) 
0(1/N*): -^ = * m A\ (25) 

^(1/iV 4 ): ^- = Pl C«S 1) + 3p 2 |«S 1) | 2 5n 2 4 1) , (26) 

with initial condition u^' (A = 0, Ai = 0, A 2 = 0, A 3 = 0) = u^p . Eqs. ^EB^^EW provide 
point symmetries of Eq. while Eq. j[2h]) is a generalized symmetry of Eq. 

Proof. The proof is done by a direct computation by taking into account the results con- 
tained in Theorem [TJ 
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Inserting Eq. (|22|) in the first symmetry (|20|) we get the following determing equations: 



0(1/N) : 



~(i) 

sin k, u\ , 



G(l/iV 2 ): ^ + 



2j9 : 

9Ai 2p 



— - (sin k ill 1 ' — iMi cos k 5 n2 u < [' > 
2rr V 



(27) 
(28) 



0(1/N 3 ) : 



<9A ' <9Ai <9A 2 
i 

2~f 



M 2 

sin k Ug — iMx cos k 5 n2 u^ H — ^- sin k S^u^ ) + 



i 



+ — ( — isin/t sm(2K)u^pu^ + Ml sin nu^'d^u^' ) + 
3i 



r,Wt „~(0) 



2p z 



• 3 |~(1)|2~(1) 

sin k\u\ \ u\ , 



(29) 



C(1/AT 4 ) 



9m 



+ 



<9w 



dA 9Ai 
i 

= 2^2 



"o 1 ^ 9m! 1 ■* SMI 1 '' 



<9A 2 



9A, 



~(1) • Ayr X ~W i ■ x2 -(i) 

sin k u\ — iMi cos k o n2 u\ H — sin k o n2 w-[ - 



iM? r3 _m . 
■— -cosK<5i< J + 



. , , /-(1)~(2) . ~(1)~(2)\ . 

-l sin k sm(2/«j lu\ 'u^ ' + U2 u\ ' J + 
+Mi sinK (^u^S n2 u^ + u^5 n2 u^ + u^S^uf^ 

-iM 2 COS KU^S^uf] + 



3i 



2 V L 



-iMi cos k sin 2 k {u^) 2 5 n2 u^ ' + 



(30) 



Let us consider Eq. (f2"7f) ; by the reparametrization A 1— > 2p 2 X/ sin k, Eq. (12T|) is equivalent 
to Eq. ( J23l) . This is the first point symmetry of the dNLS (fT2l and it corresponds to a 
phase symmetry. 

Eq. (128"]) has to be split into the following equations to avoid secularities: 



_ = _oos««y Ba «i , 



<9A 2p 



1 . ~(i) 
- sm k Uo . 



(31) 
(32) 
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From Eq. ( 1321) we see that depends on A as uf . Eq. (!3B provides the second point 
symmetry (124"]) of the dNLS f[T2|) . corresponding to translations w.r.t. the index after 
the reparametrization Ai 1— > 2p 2 Ai/(M 1 cosk). 

From Eq. f[2"9"]) . taking into account Eqs. (I27|32p and the secularity conditions, a straight- 
forward algebra and the reparametrization A2 1— > 4p 2 pi\ 2 /(Mf sin k) leads to 

r2 ~(1) , ~(1)|~(1),2 

1 ~d\ 2 ~ = Pl n2 1 P 1 ^ 1 ' ' 

which leads to Eq. ( 125]) thanks to Eq. (112p . Eq. ( 125]) means that the dNLS (fT2l) is invariant 
under translations w.r.t. the index m<i- 

Finally, Eq. (130]) gives Eq. (126]) after a long computation by taking into account Eqs. 
(I27|28ll29p . In this last case the reparametrization of the group parameter reads A3 1— > 
12p 2 p 1 X 3 /(MfcosK). 

□ 

A computation up to order 1/iV 4 , similar to the one just done for the symmetry (j2*Ul) . 
shows that the multiscale expansion of the second generalized symmetry f[2"TT) of the lpKdV 
(JH]) gives the same symmetries (I23|I24|I25||2T)]) . after suitable reparametrizations of the group 
parameters. 



5. Concluding remarks 

In this paper we have considered the multiscale expansion of the spectral problem and 
of the symmetries of the partial difference integrable lattice potential KdV equation. By a 
proper choice of the spectral problem of the lpKdV we have been able to derive from it the 
spectral problem of the reduced equation, a nonlinear Schrodinger equation. We then did 
the multiscale expansion of two generalized symmetries. A generalized symmetry provides 
us with the point and generalize symmetries of the nonlinear Schrodinger equation. At 
each order of the multiscale approximation, we get by reduction from the request that 
no secular condition exists, a higher order symmetry. The same calculation for other 
generalized symmetries do not provide anything new. All the information concerning the 
whole hierarchy of generalized symmetries for the NLS is contained in the first generalized 
symmetry for the lpKdV. 
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